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Abstract. 

We describe a realistic scheme for coupling atoms or other qiiantum emitters 
with an array of coupled optical cavities. We consider open Fabry-Perot 
microcavities coupled to the emitters. Our central innovation is to connect the 
microcavities to waveguide resonators, which are in turn evanescently coupled to 
each other on a photonic chip to form a coupled cavity chain. In this paper, we 
describe the components, their technical limitations and the factors that need to 
be determined experimentally. This provides the basis for a detailed theoretical 
analysis of two possible experiments to realize quantum squeezing and controlled 
quantum dynamics. We close with an outline of more advanced applications. 
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1. Introduction 

The interaction of light and matter is of central importance to research in quantum 
optics \T\. The strength of this interaction depends on the intensity of the light at 
the location of the atom. To enhance the interaction, photons can be trapped in a 
high-finesse cavity of small volume leading to high intensity fields inside the cavity 
even with small photon numbers. This arrangement forms the central paradigm of 
cavity quantum electrodynamics (CQED) [2], |3l 1^. Early experiments in CQED 
used atoms that where dropped through the cavity and hence passed briefly through 
the region of strong interaction [5]. With the advent of laser cooling and trapping 
techniques, it became possible to position atoms or ions at a desired location within 
the cavity and to achieve "strong" atom-photon coupling, i.e. a coupling rate faster 
than both the atomic decay rate and the cavity decay rate. Phenomena such as 
vacuum Rabi splitting [6J, photon blockade [71 E] and single photon sources [9] have 
now been observed. 

So far, experiments in CQED have almost exclusively been carried out with 
single cavities. In recent years however, it has become possible to fabricate 
microcavities in regular arrays in various settings [TUl HH [T21 [131 IH HSl [ISl HZj- Due 
to their small size, these cavities all have small mode volume and correspondingly a 
strong interaction between light and matter. Indeed the "strong coupling" regime 
has already been demonstrated for several of these devices [HI [191 [E] • 

An important next step is to find a way to couple these cavities so that they 
form an array in which photons can tunnel from one cavity to another. Recently, 
various approaches have been put forward for using such arrays of coupled cavities 
as a platform for quantum simulators [201 1211 [221 1231 [21] ■ These include a scheme 
for simulating the Bose-Hubbard Hamiltonian [201 [251 |26l |2Z] , models of interacting 
Jaynes-Cummings Hamiltonians [221123], and an effective spin Hamiltonian [23]. The 
phase diagrams of these models have been studied [291 1301 [311 132] and the existence 
of a glassy phase has been predicted [29] . 

In this work, we describe the design for a practical device, in which atoms 
interact strongly with high finesse microcavities and, at the same time, photons 
can tunnel with low loss from one cavity to the next through an interconnecting 
waveguide chip. The main components of the device are shown schematically in 
Fig. HI 

The paper is organised as follows. In section |2} we describe the design of the 
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device in detail and specify some realistic parameters. In section |3] we calculate the 
spectrum of one composite cavity. We then show (section |4]) that the dynamics of 
the device is well approximated by a Jaynes-Cummings Hamiltonian. We describe 
experiments on a small scale to demonstrate two types of quantum simulator using 
coupled cavity arrays. As a first example (section [s]) we show that a lattice of 
interacting Jaynes-Cummings Hamiltonians could be implemented in our device 
and analyse its steady states in a driven dissipative regime for a two-cavity setup 
[331 [311 ESI ES]- In a second example (section [6]), we show that our device is also 
suitable for implementing effective spin Hamiltonians and we study the dynamics in 
a two-cavity setup. We discuss briefly a variety of applications that can be based on 
these capabilities but go beyond the two cavity setting. Finally we summarize our 
results and give an outlook on future perspectives in section [8j 

2. Device 

We envisage a device with quantum emitters in several separate microcavities, 
coupled by waveguides on a coupling chip as illustrated in Fig. [T] The microcavities 
are based on the design described in [19]. Hemispherical micro- mirrors, wet- 
etched into a silicon chip [TJl HH] and reflection-coated, form one side of a set 
of Fabry-Perot microcavities. These cavities are closed by the plane reflection- 
coated ends of waveguides integrated on a photonic chip. The other end of each 
waveguide is also reflection-coated to form a second set of optical resonators. Each 
waveguide/microcavity pair, coupled together through the shared mirror, forms a 
composite cavity. Henceforth, we will sometimes simply call this a cavity, and we will 
speak of a microcavity or a waveguide cavity if we wish to distinguish the component 
parts. The microcavities are open in the transverse direction, giving access to lasers 
that trap and manipulate atoms at the position of maximum interaction with cavity 
mode. This design benefits from the intrinsic scalability of microfabrication to 
achieve controlled nearest-neighbour coupling of many optical cavities. 

In the following discussion, we assume a wavelength of A = 780 nm, that being 
the D2 line of rubidium-87 atoms whose amplitude decay rate is 7 = 27r x 3 MHz. 
In a microcavity of length Lc, the maximum coupling rate qac between a two-level 
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Figure 1. Outline of the waveguide chip. W and C refer to the waveguide cavities 
and the microcavities respectively. Circles indicate the positions of active heating 
elements. These tune the resonant frequency of each waveguide, and adjust the 
coupling between adjacent waveguides. Rectangles indicate the coupling regions. 
See Table [2] for a description of the symbols used. 



atom and the standing-wave field of one plioton is 



wliere wq is tlie 1/e^ radius of the Gaussian intensity profile at the waist and c is the 
speed of light, microcavities of the type envisaged can have mode waists down to 
2/xm, while the length of the microcavity Lc is in the range of 10 — 100 /im. Therefore 
the coupling rate is of order 27r x 0.1 — 1 GHz, and hence qac ^ 7- The reflectivity 
of each mirror is given hy Ri = 1 — (Tj + Aj). Assuming that the power transmission 
and absorption coefficients fulfil Ti,Ai ^ 1, the cavity field amplitude decays at a 
rate 

c^c ^ Tc + Ac + Tew + Acw 
Kc = —— with ~ ^ , (2) 

ZLc I 

The subscript C denotes the concave microcavity mirror, while the subscript CW 
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Definitions and Symbols 



Symbol 


Definition 


Typical values 


Lw, Lc 


length of the cavities 


Lc = (10 - 100) //m 
Lw > (10 - 20) mm 


n 


amplitude reflection 
coefficients of the mirrors 




Ri = Inp 


Intensity reflection 


99.95% > Rc = 99.9% > 99.0% 




coefficients of the mirrors 


99.95% > Rw = 99.0% 
Rcw ~ 98% 


FSRi = no/ Li 


free spectral range 


FSRc 27r X (0.2 - 2) THz 
FSB^xr 27r X fl — 2"! GHz 


Fn- 

i-VRcRcw 

TP TT 


cavity finesse 

VV Ci V VyCi UllV^Vy XiXlvyOOvy 

resonance linewidth 




l-^JRwRcw 
Ki = FSRi/2Fi 






microcavitv mode waist 


itir' = 3 — 5 /ym 

kAJ (J \J '>_/ LAi A i J. 


WW 


waveguide mode size 


ideally matched to wc 


„ _ / 3cA27 

Qac - y T,2^iLc 


atom-cavity coupling 


27r X (0.1 - 1) GHz 


7 


atom amplitude decay rate 


Rb: 7 = 27r X 3 MHz 


A 


wavelength 


Rb: A = 780 nm 


Jww 


tunnelling rate between ad- 
jacent waveguide resonators 


< Jww < 27r X 2 GHz 



Table 1. Definitions of relevant symbols. All frequencies are angular: e.g. 
27r X IGHz means 6.28 x lO^s'^ 



denotes the plane coupling mirror between the cavity and the waveguide. As we 
intend to fabricate the concave mirror by isotropic etching of silicon, we expect 
the losses due to surface roughness to be of order Aq ~ 10^"^ without additional 
polishing. The coupling mirror, on the other hand, can have losses on the order of 
Acw ~ 5 X 10^^, assuming the waveguide facet is super-polished and a dielectric 
mirror formed by ion-assisted deposition is used. The decay rate of the microcavity 
could then be made as small as ~ 27r x O.OlGHz for a 100 //m-long cavity. This 
length will give enough space to permit external optical access to the atoms. However, 
in order to couple effectively to the waveguide resonator, it is desirable to increase the 
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transmission of the coupling mirror so that Tew ^ Tc + Ac- Therefore in practice, 
Kc > 27T X 0.1 GHz. This does not necessarily imply higher loss for the composite 
cavity since photons that go through this mirror enter the waveguide and are not 
necessarily lost. 

Similarly, the field in the waveguide resonator decays at a rate 

'^^w , ^ + Aw + Tew + ^cw + Kw 
i^w = wi— with ^w ~ ^ , (3) 

where Lw is the optical length of the waveguide cavity (refractive index times 
physical length) and Kw accounts for the waveguide propagation loss over one 
round-trip. While there is little use of integrated waveguides at 780 nm, we find 
in the literature [2Z1 that propagation losses can be less than 0.02dB/cm for 
wavelengths around 800 nm in polymer waveguides. This platform offers all the 
necessary technological components for our device. We estimate that to achieve the 
necessary coupling lengths and separations between the waveguides without incurring 
additional bend losses, the waveguide resonator will have a length of Lw > 1 cm. We 
will therefore conservatively consider Lw = 2 cm, which gives a fractional round-trip 
loss of 1.8 % and a corresponding decay rate of kw > 27r x lOMHz. 

The coupling from the microcavity into the waveguide resonator depends on the 
transmission Tew, but also on the spatial overlap of the waveguide and microcavity 
modes, which will need to be optimised experimentally. In principle, this can reach 
unity, and we estimate that at least 90% should be achievable in practice. The 
photon tunnelling rate between microcavity and waveguide can easily exceed the 
free spectral range of the waveguide cavity. In this case, it is appropriate to consider 
the eigenmodes of the composite cavity, rather than viewing the microcavity and 
waveguide cavity as individual devices. We calculate the spectrum of one such 
coupled cavity in section |3| 

Each waveguide is coupled to its two nearest neighbours by short regions of 
evanescent field overlap to produce a tunnelling rate Jww- The maximum value 
Jww"^ ^ determined by the maximum field overlap and the length of the coupling 
region, can be as high as 10^°s^^. This overlap can be tuned using thermal phase 
shifters, one near each end, to move the standing wave without changing the cavity 
length. In this way, Jww can be tuned in principle over the range — J^^''. 
In practice, a small amount of crosstalk will be unavoidable because of scattering, 
tuning noise and linewidth effects, but we estimate that j'^^^ < 10~^ x J^^K The 
maximum coupling strength must therefore be chosen judiciously to ensure access to 
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Figure 2. Calculated reflection spectrum of two composite cavities having 
evanescent-wave coupling between their waveguides. A pair of thermo-optics phase 
shifters is adjusted to give a phase shift 50 between the standing waves in the two 
waveguides. The system is pumped on waveguide number 1. At 6(j) = 7r/2, the 
overlap integral of the standing waves is zero, so the coupling is switched off. 
Consequently, there is no normal mode splitting and no power is seen in waveguide 
number 2. When the phase shift changes to Scf) — or tt, the coupling is maximised 
and therefore so is the mode splitting. This calculation is based on an expanded 
version of the transfer matrix method described in |47j . 



the desired range of coupling. The thermo-optic phase shifters also allow the optical 
length of each waveguide cavity to be tuned at a rate of up to a few wavelengths per 
ms. Figure [2] illustrates how two coupled cavities exhibit a normal mode splitting 
and shows how the evanescent-wave coupling can be adjusted by shifting the phase 
of one standing wave relative to the other. In particular, the coupling is switched off 
for a relative phase of 7r/2. 

3. The composite cavity spectrum 

As our starting point, we consider the spectrum of a single composite cavity. Several 
such cavities are to be concatenated in our setup to form the coupled-cavity array. 

Suppose that we couple light into one of the cavities, through the top waveguide 
mirror as drawn in Fig. [T} Without the curved microcavity mirror at the bottom. 



Arrays of waveguide- coupled optical cavities that interact strongly with atoms 



8 



the field refiected by the waveguide cavity is given by 

Er,W = --. ^-7— , (4) 

where E^ is the incident field, rw and rcw are the amplitude reflection coefficients 
of the two mirrors, and 20vk = 47rLvi//A is the round-trip propagation phase. Here 
we have assumed lossless mirrors for simplicity. Resonance occurs when (pw = rmr 
with m being an integer. Similarly, a travelling field -Ejn.c that is incident on the 
microcavity produces a refiected field in the waveguide given by 

Er,c ^ raw - rce'^''- ^ e 

where 20(7 = AttLc/X- This ratio defines the refiection coefficient rc and phase shift 
6 at the bottom end of the waveguide cavity when the microcavity mirror is in place: 

~2 ^ r^w + i^c- '^rcwrc cos(20c') 
^ 1 + r^^rl - 2rcwrc cos(20c') ' 

e = arctan ( {rlw ' ^VcM'^M \ 

Vcwirl + 1) - rcirlw + 1) cos(20c) / ' 
On replacing rcw in Eq-Q by rc^^^, we obtain the refiected field from a composite 
cavity: 

„ rw - rce'^'^^'t'^'^ 

which is in resonance when 

+ 9 = 27m, (8) 

for any integer n. This set of equations describes the system fully. It cannot generally 
be solved analytically, but we can determine when the refiected intensity becomes 
zero. This is of interest because the spectrum is generally sensitive to the presence 
of atoms under this condition. Zero refiected field is achieved when the microcavity 
can be tuned to make = r^. Such tuning is possible when the following inequality 
is satisfied: 

/ \ 2 2 2 

|7T.mm|2 VC — T^CW) ^ 2 ^\rpmax\2 '^C^''^CW ((^\ 

' = (l- rcrc,.)^ - - ""-^ I = 1 + (rcrcwY ' 
and occurs when the cavity round-trip phase satisfies 

2^c(opt) = arcco= 2rc„,rc(rJ, - 1) j ' ^'"^ 
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Figure 3. Behaviour of composite cavity with cavity lengths Lc — 156 /im, 
L\Y = 15.6 cm, and power reflectivities Rc = 0.999, Rcw = 0.98, and Rc — 0.99. 
Top left: Calculated reflection spectrum of the composite cavity versus detuning 
of the laser. In the absence of coupling, each separate cavity has a resonance at 
Auj = 0. Here they are coupled to produce normal modes split by 2.5 x lO^^s"^. 
Dashed curve: f^. Top right: reflected intensity near one of the central resonances. 
The laser frequency is fixed at Aw = 0, while the (uncoupled) resonance frequencies 
of the two constituent cavities are varied. The green circle indicates the point of full 



contrast where the reflected intensity is zero, as calculated from Eq. 10 Bottom: 
Circulating field strength inside the waveguide resonator (left) and the microcavity 
(right), as a function of their detuning from the laser. Green circles as as above. 



This determines the phase shift 6, and the resonance of the whole cavity is then 
ensured by adjusting the length of the waveguide to satisfy Eq. 

The semi-logarithmic graph on the top left of Fig. [3] shows how the reflected 
intensity varies with a change Au in the angular frequency of the laser for a cavity 
with L]Y = 15.6 mm and Lq = 0.156/im, and with mirror reflectivities of Rc = 0.999, 
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Rcw = 0.98, Rw = 0.99. Without any couphng, both constituent cavities are 
resonant at Au = 0. In the reflection spectrum we see a series of dips, symmetrically 
disposed around Au = 0. The two central dips represent normal modes of the 
composite cavity. These would be degenerate in the absence of coupling but are 
split apart here. Further away from this doublet, the resonances are not so strongly 
coupled to the microcavity and their spacing approaches the free spectral range of the 
bare waveguide cavity. The horizontal line at 0.99 shows the reflection coefficient Vyy 
of the input mirror, while the dashed line indicates the value of f^. These are most 
nearly equal at the second reflection dip on each side, making those dips the strongest 
ones. The width of each resonance can be understood in a simple way because the 
waveguide resonator is very much longer than the microcavity. Consequently rc and 
6 are essentially constant over any given resonance line of the coupled system. By 
analogy with Eq. ^ this gives the cavity damping rate as 

k ^ 77^^ with ^ ^ 1 - rwfc, (11) 

The top right graph in Fig. [3] explores how the reflected intensity near one of the 
central resonances depends on the detuning of the two constituent cavities. We see 
that the reflection goes through a minimum whose position, indicated by the small 



circle, is as expected from Eq.(lO). The two lower graphs show how the circulating 
power inside the waveguide resonator (left) and the microcavity (right) varies with 
the same detuning of the constituent cavities. In particular, we see that the condition 
of minimum reflected power (again indicated by a circle) corresponds closely with 
having maximum power inside the microcavity, as required for good atom-cavity 
coupling. At this point, the energy density in the microcavity is ten times that in 
the waveguide (though the total energy is ten times less because of the hundredfold 
disparity in lengths). 



4. Approximation by a Jaynes-Cummings lattice Hamiltonian 

In this sections we describe a series of approximations that allow us to obtain a 
description of our system in terms of atom photon interactions of Jaynes-Cummings 
form in each cavity together with tunnelling of photons between adjacent cavities. 
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4.1. Composite cavity Hamiltonian 

The Hamiltonian for the field in one composite cavity is 

a 

where the index j labels the particular cavity and index a labels its eigenmodes. The 
frequencies cUq are those of the eigenmodes discussed in Sec. [3] 

4-2. Photon tunneling between adjacent waveguides 

The rate Jww for photons to tunnel between two adjacent waveguides can be tuned 
over a wide range. Let us choose to make it small compared with the waveguide 
free spectral range, which also ensures that |Jvi^vk| ^ \oJa — 0Ua'\- In this regime, 
we can write the coupling as a tunnelling term between resonant modes of the two 
composite cavities, 

a,/3 a 

When the cavity network only contains photons that are near-resonant with the 
normal modes a^, and when the Rabi frequencies of atom-photon coupling are also 
^ \ua —oJa'l, only these modes are populated and Hww can be simplified further to 
read, 

Hww ~ —Jao,ao {^\oio'^j+l,ao + H.C. j . (13) 

The tunnelling rate Jao^ao is related to Jww- Yet, since Jww is harder to determine 
experimentally than Jao,ao) we do not specify the conversion here. 

4-3. Atom-photon coupling 

The circulating fields in the microcavity and waveguide, Ec and Ew respectively, 
are given by, 

E =E tcwtwe''^^ 

*^ 1 - rcTcwe'^''^'^ + rcrwd'^''''*"^^'^^'' - rcwrwe'^^*'^ 
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where Ei^ is the input field from the top end of the waveguide, see figure [TJ Since one 
photon in the composite cavity has a total energy content of hw, the absolute values 
of the circulating field amplitudes per photon in the microcavity and waveguide are 



^ ^ 2eo7rw^y LclEcf + Lw\Ewf ^^^^ 



^1 2eo7iwl\] Lc\Ecf + Lw\Ewf' ^^^^ 
Hence the atom-photon coupling at an antinode of the microcavity is 

(18) 



9 = 9 AC 



\Ec\^ 



^ \Ec\' + ^\Ew\' 

where qac is given in eq. Q for a 2-level atom. This formula for a two-level atom 
applies because we are considering the closed cycling transition 5^S'i/2(-F = 2,mp = 
±2) ^ 52P3/2(F' = 3,m'p = ±3) on the D2 line of ^^Rb. 

4-4- Jaynes-Cummings lattice 

With the above approximations, the dynamics of atoms coupled to the composite 
cavity normal modes be described by a Jaynes-Cummings lattice model. 

Dropping the index on the field operators, the Hamiltonian for composite 
cavities reads, 

N N 



-f^JCarray = ^ CTj^CT- + UJc^ a',a 



3=1 3=1 



N 



i=i 

7V-1 

Here, uc = cuao, o"~ = \9j){ej\ is the transition operator between the excited state 
\ej) and ground state \gj) of the atom in cavity j and is the annihilation operator 
for photons in mode cto of that cavity. The losses for the system described by this 
Hamiltonian arise through spontaneous emission from the excited states of the atoms 



at rate 7 and photon loss from the normal modes a,- at rate k given by Eq. (11). 
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5. Spectroscopy for a driven Jaynes-Cummings array 

We now consider an array of composite cavities connected to each other by nearest- 
neighbour coupling and excited at the end of one waveguide by a laser tuned to 
the normal mode ai. The resonant pumping is described by an additional term 
Hd = + ill the Hamiltonian whereas photon leakage and spontaneous 
emission losses are taken into account by Markovian damping terms. The dynamics 
of the resulting driven dissipative system is then described by the master equation, 

3 

+ k ^^(2ajpaj — aja^p — pa^aj), (20) 
j 

where H = Hd + -ffjcarray For one cavity, a driving amplitude rj gives a steady 
state energy in the cavity of hu){a^a) = hu)\ri\'^/ (2k^). Also, the ratio of energy in the 
cavity to input power P is 2 {\Ec\'^Lc + \E]y\'^L]y) /(\Ei^\'^c), where Ec and Eyy are 
the circulating fields in the cavity Hence, the power P of the driving laser is related 
to the driving amplitude t] by, 

Ir^P E^c 
^ = ^'^2i\Ecncl\E^\-L^y 
Below, we calculate the spectrum and photon statistics that can be observed at 
the output ports of the coupled waveguides [33], [Ml EH |36]. These properties 
correspond to the most straightforward experiments that might be made using such 
a cavity array As we shall see, they nonetheless reveal interesting physics including 
significantly entangled photon output states. 



5.1. Composite cavity modes with strong coupling and high single atom cooperativity 

Since we are interested in having large atom-photon coupling and small photon losses, 
both from the cavities and from spontaneous emission, we wish to work with a mode 
that has high single atom cooperativity, Ci = g"^/ {2k'-)) and is resonant with the 
atomic transition. The required combinations of microcavity and waveguide lengths 
are found by numerical optimisation. In Fig. |4] the cavity resonances are indicated by 
vertical lines. The resonance near (5^ = 0, with the highest cooperativity (solid line) 
of Ci = 17.3, coincides with the atomic transition frequency. The cavity parameters 
{Lc = 156.05 /im, Lw = 20.000 mm, Rc = 99.9%, Rcw = 98.0% and Rw = 99.8%) 
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Figure 4. Composite cavity properties versus laser detuning. The atomic 
transition frequency is at buj = 0. The cavity parameters (Lq — 156.05 /xm, 
Lw = 20.000 mm, Rc = 99.9%, Rcw = 98.0% and Rw = 99.8%) are close to 
those of Fig.[3j Vertical lines: frequencies of reflection minima due to composite 
cavity resonances. Solid line: single atom cooperativity C\ = (7^/(2k7). At its 
peak value C\ = 17.3. Dashed line: coupling relative to atomic decay rate, (7/7. 
Dotted line: coupling relative to cavity decay rate, gj^. 



are very close to the parameters used in Fig.[3j This resonance corresponds to the 
strong reflection dip in Fig. [3] on the high frequency side of the central doublet. The 
atom-photon coupling g and photon loss rate k are g = 0.3304x g^c = 27rx 33.04 MHz 
and K = 27r X 10.6 MHz. 

5.2. Steady-state of two- site Jaynes-Cummings array 

Figure [5] shows the steady-state behaviour when two of these cavities are coupled 
together at a rate J = 0.2 x g^^c = 2tt x 20 MHz. Each contains a rubidium atom 
and cavity number 1 is driven by a laser with 77 = Q.lxgAc = 27rx 10 MHz. The atoms 
are resonant with the normal modes aj. Figure [5^ shows the number of photons in 
each cavity: ni = {a\ai) and n2 = {a\a2). The four resonance lines in the spectrum 
correspond to the four singly excited eigenstates of the two-site Jaynes-Cummings 



Hamiltonian, i.e. Eq.(19) with N = 2. Here, the resonances at 5a; ~ ±27r x 44 MHz 



correspond to states where the excitation is more likely to be found in one of the 
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cavity modes whereas for the resonances at ~ ± 27r x 24 MHz the excitation is 
more hkely to be in the atom than in the cavity field. 

In figure |5}d we study the photon density correlations g^^\l, 1) = {a\a\aiai) / nf, 
g^'^\l,2) = {a\a\a2ai) / nin2 and g^'^\2,2) = {a\ala2a2) / n^. On the two outer 
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Figure 5. Steady state of a driven array of two Jaynes-Cummings cavites for 
g = 2TTX 33.04 MHz, r; = 27r x 10 MHz, J = 27r x 20 MHz, k = 27r x 10.6 MHz and 
7 = 27r X 3MHz. a) shows ni and na, b) shows g'^'^'>{l,l), g^^^{2,2) and g'^'^'>{l,2) 
and c) shows the entanglement between modes ai and 02 as quantified by the 
logarithmic negativity E'at. 
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resonances {6u) ± 27r x 44 MHz) we see small dips below unity, indicating photon 
anti-bunching in both cavities, g'^'^\l,l) < 1 and g^'^'>{2,2) < 1, as well as anti- 
correlations between photons in distinct cavities, |5f*^^^(l, 2) | < 1 . This indicates 
that the coupling in both cavities is strong enough to generate an optical nonlinearity 
that converts coherent classical input light into a manifestly non-classical state of 
the cavity photons. For the resonances at 6uj ^ ±2tt x 44 MHz the state of the 
two cavity modes is approximately a superposition of one photon in cavity 1 with 
cavity 2 empty and one photon in cavity 2 with cavity 1 empty which gives rise to 
the anti-correlations shown by g^^\l, 1), g^'^\2, 2) and ^''■^■'(l, 2). Since cavity 2 is not 
directly driven, it experiences a lower intensity of incoming photons than cavity 1 
and hence exhibits slightly stronger anti-bunching for the same nonlinearity. For this 
state of the cavity modes, the photons in the two cavities are expected to become 
entangled. This is confirmed by figure [S}^ which shows the entanglement between 
modes ai and 02 as quantified by the logarithmic negativity [39] . 

En = log, (tt^I^?! , (22) 

where ppt is the partial transpose of the reduced density matrix of the two modes ai 
and a,. 

On the inner resonances {600 ^ ± 27r x 24 MHz), by contrast, only photons 
involving both cavities are anti-correlated whereas photons within each cavity bunch. 
This indicates that photons prefer to stick together in either of the cavities and the 
associated state contains superpositions of two or more photons in cavity 1 with 
cavity 2 empty and the reverse configuration. Hence photons in the output of one 
cavity tend to come in bunches whereas if one cavity emits a photon the other 
cavity is unlikely to emit at the same time. Consequently we also find entanglement 
between the photon modes for these resonances although less than for the resonances 
at 5a; ~ ± 27r X 44 MHz which feature higher photon densities. 

Close to 6u = 0, the photons in the cavities show pronounced bunching. This 
emerges due to the nonlinear spectrum of our device which here causes the laser 
drive to be detuned with respect to single photon transitions but resonant with 
multi photon transitions. In practice this bunching will however be hard to observe 
since the photon densities generated in these multi photon processes are vanishingly 
small as apparent from figure [5^. 

The findings shown in figures |5]d and c clearly demonstrate the non-classical 
nature of the light fields generated in our device. 
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6. Effective Spin Hamiltonians 

In this section we describe a second experiment that could be performed using such 
a device. We show how effective spin-spin interactions, as proposed in [2B|, can be 
implemented. Here, each cavity interacts with one Rubidium atom, whose energy 
levels a, b, e associated with the D2 line form the lambda structure depicted in figure 
[6j These levels are coupled both by the cavity photons (with couplings ga and gb) 
and by external laser fields (with angular frequencies z/^, Uh and Rabi frequencies fla, 
Qb)- Importantly, the transition \a) o \b) (5^5'i/2(-F = l,mF = 1) o 5^5'i/2(-F = 
2,mF = 1)) is dipole forbidden and level \b) is thus metastable. Under conditions 
spelled out below, the dynamics can be constrained to the subspace formed by levels 
\aj) and \bj) of each atom and we can identify \aj) with spin down, = \aj), and 
\bj) with spin up, = \bj). 




OJb 



a) 



Figure 6. Lambda level structure for the D2 line of *^Rb atoms. The state \a) 

represents 5^S'i/2(F = l,mp = 1), |6) stands for 5^S'i/2(i^ — 2, nip = 1) and |e) is 
5^^3/2(^' = 2,mF' — 2). The transitions \a) -H> |e) and \b) O |e) are driven both 
by the cavity field (green arrows) with coupling strengths ga and gb, and by external 
coherent fields (blue arrows) with Rabi frequencies Qa and Of,. Sa = uJe — and 
Sb = u!e — {ujb — 5) — ujc- State \b) is stable as there is no electric dipole moment 
for the transition \a) <-> \b) 
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6.1. Outline of the approach 

Following the arguments presented in section [3} we assume once again that the atom 
in the j*'* composite cavity only couples to one normal mode aj (again, we skip the 
index ao), for which we maximise the single atom cooperativity as in section [s] The 
photons tunnel between adjacent composite cavities at a rate J. The Hamiltonian 
for this system can thus be written as 

(23) 



H 



Ha + He + Hac where 



N 



Ha =Y^ Uelcj) {cj I + ujb\bj) {bj 

N N 

He = coc ^ (ijO-j — J ^ (ojaj+i + ajaij_^_^^ and 

N 

j=l x=a,b '- ^ 

Here, Ha describes the atoms. He the cavity field and Hac the interaction between 
atoms, lasers and cavity field. Also, Ue is the a — e transition frequency and Ub the 
a — b transition frequency. 

The Hamiltonian Hq can be decomposed into non-interacting collective photon 
modes. He = J2k ^kO'l^^k, where Uk = ujc ~ 2 J cos k and 



+ 9x0.] \ej){xj \ + H.c. 



(24) 
(25) 
(26) 



EliSin(/ci)c 



N+i 



with k = -j^^ and / = 1, 2, . . . A^. It is helpful to move to an interaction picture, 

H{t) = e'^''\H - /Jo)e"'^°* (27) 
where H is given in equation (p3| and Hq reads. 



N 



Ho = Yl i^e\ej){ej\ + {ub - S)\bj){bj\) + ^ujkalak, 
j=i k 



(28) 



with the value of 5 to be chosen later. In this picture, the Hamiltonian (23) becomes, 

N 

H{t)=5j2\bj){bj\+ (29) 



N 

j=l x=a,b k 



"e^^^* + ^,,,-fce^''^^*afc) |e,)(x,-|+H.c. 



i5?xt 
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6) — z/fe, as illustrated in Fig.[6| 



Up 



where Aa = Ue - l^a, A^ = Ue - {Ub 

and 6l = LOe — (wb — 5) — ook- The coupling constants ga,j,k and gbj,k are related to 
the couplings ga respectively gi, via gx,j,k = a / t^ttt sinfA;?)^^, for x = a,h. 



N+ 



A judicious choice of rotating frame is 5 



\{ya 



which ensures that 



"k 



Ab = 6l — Aa for all k and allows a convenient separation of fast rotating terms 
from near-resonant terms in what follows. In this rotating frame, a second order 



adiabatic elimination of the excited levels \ej) and the photons aj, see Appendix A 
yields the effective spin-^ Hamiltonian, 



N 



spin 



where the effective transverse field B reads. 



(30) 



B 



+ 



5_ _ 
2 ~ 

I a 



8Al 



A, 



SAl 



A„ - 



I a 



2Ab 4(A„ - Ab 



E 



\9b,j,k\ 



A, 



+ 



\9a,j,k\ 



61 - Ab 



2Aa A{Ab - A, 



k 



in 

\ya,j,k I 



+ 



\9b,j,k\ 



6t 



Aa 



and the coupling constants are 

^a^h Sr^ 9b,j,k9a,l,k 



4A„A, 



E A'' _ A 



K 



\2 Z-^ 



4A 



9b,i,k9b,j,k 
61 



+ 



w 

4A? 



for X, y 



a, b, 



k-^a 
O 

1 for X = a,b and 



9a,j,k9a,l,k 



The Hamiltonian (30) is an appropriate description provided that lA^; — ^ FSRw 

l2A^ 



<^ 1 ioT x,y = a, 6, see 



conditions (A.5). Furthermore, all the eigenmodes of the composite cavity should be 



sufficiently detuned from the |e) — ?■ \a) and |e) — )■ |6) transitions to avoid Purcell- 
enhanced atomic relaxation. 

The sums over all photon modes, depend on the number of cavities in the 
array. As a starting point for experiments, we focus here on the N = 2 case. The 
explicit expressions for this case are given in Appendix B 
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6.2. Effective spin dynamics for two coupled cavities 



As an example of dynamical evolution under this effective spin Hamiltonian, we 



consider once again the two coupled cavities described in Sec. 5.2, each containing 



a ^^Rb atom. The atom in cavity number 1 is prepared in state while the 
other atom is placed in state [02) to form an initial effective spin state |ti7 4'2)- 



The subsequent evolution of this state under the Hamiltonian of equation (30) is 



illustrated by the dotted and dash-dotted lines in Fig. [7} which plot the probabilities 
of the spin-up and spin-down states in cavity 1: P^i = Tr(|ti)(ti| p) and 
Pl^i = Tt{\Ii){Ii\ p), p being the state of the effective spin system. We see that 
the spin oscillates in this case with a period of approximately 70 /xs, determined by 
the coupling constant K in equation (30). The symmetry of the problem ensures 
that P^^2 = Pi,i and Pi^2 = P\,i- The values of ga and gb are derived from equation 



(18) (via equations ([TJ, ( [m] ) and (15)), together with the appropriate weightings 
(l/\/2 and l/\/6) relative to the cychng transition. The values of these and the 




15 

time (|as 



Figure 7. Evolution of effective spin using the two-atom, two-cavity system 
described in Sec. |5.2[ The cavity details are given in the caption of Fig. |4] Other 
parameters are J = 27r x 100.0 MHz, g„ = 27r x 23.36 MHz, gb ^ 2tt x 13.49 MHz, 
k^2n X 10.59 MHz, = 27r x 166.67 MHz, = 2n x 394.16 MHz, Aa ^ 2tt x 
5.000 GHz, Afc = 27r X 11.825 GHz, 6a = 2tt x 11.740 GHz, (5f, = 27r x 4.915 GHz and 
(5 = 27r X 10.0 MHz. These values lead to B = 27r x 4.05 MHz, J,j = 27r x 6.188 kHz 
and Kj I = 2-k x 7.145 kHz. Dotted and dash-dotted lines show the probabilities for 



states 14,1), Iti), according to Eq.(30). Dashed and solid lines show probabilities 



for states |ai) and |6i) according to the Master equation given in Appendix G The 
damping has a significant effect. 
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other relevant constants are listed in the caption of Fig.[7| 

In deriving equation (30) we adiabatically eliminated the excited state |e) 
and therefore ignored the spontaneous emission from the atoms. In reality, this 
emission dephases the effective spin and causes loss of probability from the three- 
level system {a,b,e}. In addition, equation (30) ignores the loss of photons from 
the cavity. If a real experiment is to simulate the dynamics of a spin chain, as 




2 



q\ , , , , , , 1 

5 10 15 20 25 30 

time in \is 

Figure 8. Dynamics for effective spin system approach for Rc = 99.99%, Rcw = 
98.0%, Rw = 99.9%, Lc = 103. ^m, Lw = 20.000mm, J = 2n x 100.0 MHz, .g^ = 
2tt X 29.64 MHz, gt = 27r x 17.11MHz, k ^ 2tt x 3.10 MHz, f^^ = 27r x 166.67 MHz, 
= 27r X 394.16 MHz, = 27r x 5.000 GHz, Af, = 27r x 11.825 GHz, Sa = 
2tt X 11.741 GHz, 5^ = 27r x 4.916 GHz and S ^2ttx 10.0 MHz. These values lead to 
B = 27r X 4.05 MHz, J,_, = 27r x 9.40 kHz and Kj^i ^ 2t: x 10.85 kHz. a) Occupation 
probabilities i for the states |ti), \ai) and \bi). b) Occupation probability 
of the excited state Pe,i and expectation value for the number of photons ni. The 
initial state is |6i,a2) respectively |ti,i2). 



Arrays of waveguide- coupled optical cavities that interact strongly with atoms 22 



given by the Hamiltonian in equation (30), these rates of dissipation must be small 
enough. To test whether this is the case, we have calculated the full dynamics of the 
relevant atomic levels, including spontaneous emission and photon leakage using the 
master equation given in Appendix C The solid and dashed lines in Fig. [7] show the 
probabilities for the atom in cavity 1 to be in states a and b: Pa^ = Tr(|ai)(ai|p) 
and Pb^ = Tr(|6i)(6i|p). These should correspond to the spin-up and -down states 
of the equivalent spin model, but we see that they do not because the coherence is 
damped, leaving an incoherent mixture of states \a) and \b). 

The situation can be rectified, as illustrated in Figure|8]^a). Here, the reflectivity 
of the microcavity mirror has been increased from 99.9% to 99.99% (together with 
a slight shortening of the microcavity and a minor improvement in the waveguide 
mirror). This improvement in the microcavity mirror allows the real system to 
provide a reasonable approximation to the ideal spin evolution, albeit with some 
residual damping. Figurelstb) shows the shows the excited state population Pe,i 



on site number 1, which determines the lifetime (Pg ^ | ^ 27)^^ ~ 260 ps for 
atoms to be lost from the lambda system by spontaneous decay to the state 
5'^S{F = 2, nip = 2). Also shown is the small population ni of cavity photons. 
Significant further improvement in the coherence is possible in principle if the 
dissipation in the waveguide cavity is also reduced. In practice, this will require 
an advance in integrated waveguide technology to achieve smaller propagation losses 
than the current state of the art. 



7. More advanced experiments based on fundamental capabilities 

The experiments that we have studied in some detail in the previous subsection serve 
to demonstrate both the ability to couple separate cavities and to create effective 
interactions between the atoms that are held in each cavity. These are essential 
capabilities that provide the basis for a wide variety of more challenging experiments 
for which we would like to briefly outline some examples here. 



7.1. Preparation of entangled states by propagation 

A stringent test of the coherence properties of transport in extended chains of 
cavities is provided by the transport of quantum entanglement through the chain. 
Here one might envisage the initial creation of entanglement by controlled quantum 
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dynamics either between directly neighbouring cavities to generate entangled two- 
mode squeezed states between two cavity modes or between atomic qubits held 
in different resonators. Most promising in this context though appears to be the 
generation of entanglement between an atomic qubit and the photonic degree of 
freedom of the cavity it is residing in. Once created, the coherent coupling between 
constituents of the cavity array will lead to transport of the state of the photonic 
degree of freedom as described in a harmonic chain [32]. At a suitably chosen time 
when we expect maximal entanglement between the atomic qubit and the photonic 
degree of freedom of a distant cavity we would then probe for the entanglement either 
by quantum state tomography or by measurement of restricted sets of observables 
as described in [H]. 

7.2. Probing coherence by transport measurements 

Another approach to determining coherence in cavity arrays exploits the quantum 
coherent effect of dynamical localization [l3t HI] in a harmonic chain jl5] governed 
by a Hamiltonian of the type 

N N-1 

H = Y^ h{no + ^1 cosut)ka+cr- + J] hc{(T+(T-^-^ + (T^(r+^^), (31) 

k=l k=l 

that is, a chain of sites subject to a hopping interaction with strength c and a static 
and a time-varying field that affect their on-site energies energies via Qq and Qi 
respectively. This may be realised by a spin model or, alternatively, by a set of 
harmonic oscillators replacing cr+ (a^) by a'^ (a) respectively. In a spin model the 
variation of the on-site energies could be obtained by a time dependent magnetic 
field gradient, which induces a shift of the energy levels. For harmonic oscillator 
systems the cavity resonance frequencies would have to be shifted for example via 
modulations of the resonator lengths or of the refractive index inside a resonator due 
to the presence of a detuned two level system. 

The existence of the dynamical localization and its coherent character may be 
seen by moving to an interaction picture defined by lipii)) = e~^^^^^^k^°'k'^k \ip{t)) 
with A{t) = —Qot — {Qi/u) sin{ut) in which the time-dependent on-site energies 
vanish at the expense of introducing a time-dependent coupling strengths between 
the neighbouring sites |44j. For small coupling strength c the hopping dynamics is 
negligible over the interval [— vr/cj, vr/cj] and then we find the effective time-averaged 
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Hamiltonian 

N-l 

Hj = Y^ hcJn,/c. i^i/uj) {(^t(^k+i + (^k^k+i) ■ (32) 

k=l 

li Qq/u = n and if Vli/uj coincides with a zero of the Bessel function J„, then the 
evolution of a wave-packet becomes periodic in time, the spreading of the wave-packet 
is suppressed, and hence transport is suppressed. This expresses itself as resonances 
in the ability of the chain to transport excitation, which one can observe directly 
without the need for process tomography. In the averaging that gives rise to the 
Bessel function, it is the destructive interference of amphtudes that suppresses the 
transport. In a classical system, this interference is absent and there is no localisation. 
The strength of dynamic localization can in fact be used to infer the level of quantum 
coherence in the system without needing detailed measurements or tomography, but 
depending purely on the quality of the transport . 

7.3. Long distance entanglement 

Systems with finite correlation length, such as the 1-D Heisenberg and XX models, 
allow sizeable ground-state end-to-end entanglement, independent of the size of 
the system, provided that simple patterns of site-dependent couplings are selected. 
This phenomenon has been termed Long-Distance Entanglement (LDE) and can be 
implemented in cavity arrays |1HI 09]. The realisation of a model with spin-spin 
interactions of XX type requires local control, as described in section 6, to create 
position-dependent coupling between neighbouring sites. In particular it will be 
necessary to couple the first and last sites of the chain to the bulk of the chain 
weakly in order to obtain strong entanglement between those two sites in the ground 
state of the chain. A chain of 4 cavities suffices to demonstrate this effect. 

7.4. Adiabatic quantum state preparation 

The generation of complex quantum states may be achieved by a sequence of 
elementary quantum operations correlating pairs of sub-systems. However, under 
realistic conditions the resulting fidelity will tend to scale very badly with increasing 
number of particles and complexity of the quantum state. It may therefore be 
preferable to take advantage of the natural dynamics. For example, we might first 
prepare the system in a product state without interactions between sites, then slowly 
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switch on the couphng between neighbouring sites. As a consequence of the adiabatic 
theorem, the system will remain in its instantaneous ground state, which thanks to 
the interactions becomes a complex entangled state. This highly entangled state of 
the atomic degrees of freedom could then serve as a deterministic source of highly 
entangled photons. Indeed, the atomic degrees of freedom may be mapped onto the 
photonic mode at any desired time and photons allowed to escape into fibres. 

7.5. Non-adiabatic creation of entanglement 

The generation of entanglement between distant sites in the cavity array is 
complicated by the fact that interactions are between nearest neighbours. Distant 
particles may be entangled through a sequence of nearest-neighbour interactions, 
but the efficiency of this preparation decreases very rapidly with distance. An 
alternative is to use the natural dynamics of states that are not eigenstates of the 
system Hamiltonian. Thanks to the availability of local control in a cavity array, it 
is possible for example to prepare the atoms inside the cavities in a product state. 
For most spin Hamiltonians, such as those generated in sec 6, this initial state will 
evolve in time to become highly entangled. In particular, significant two-particle 
entanglement will build up between sites whose distance is proportional to twice the 
speed of sound multiplied with the waiting time . A basic demonstration of such 
an experiment is already possible with three sites but may be scaled to larger arrays 
as the strength of the entanglement decreases only slowly with distance. 

8. Summary and Outlook 

In summary, we have presented a practical way to realise an array of coupled cavities 
that could be used for quantum simulation. The device consists of open Fabry-Perot 
microcavities that are coupled via a waveguide chip. We have demonstrated that, 
under suitable conditions, the dynamical evolution is well approximated by a Jaynes- 
Cummings lattice Hamiltonian and that it is suitable for implementing effective spin 
Hamiltonians. We have outlined experiments, both basic and more advanced, that 
could be carried out in such a device to achieve controlled quantum dynamics. 
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Appendix A. Derivation of effective spin interactions 



The Schrodinger equation containing the Hamitonian H as in equation (29) reads, 
'^\<iI,t) = ~zH{t)\^,t). 



dr 



This equation can formally be integrated to yield 

i-t+T 

\m,t + T) = \m,t) -I dsH{s)\^,s). 
Iterating the right hand side results in. 



|^,t + T) = + 

ft+T 



n=l 



dtnH{tr, 



dtn_lH{tn- 



ti 



(A.l) 

(A.2) 

(A.3) 
dtoH{to)\^,t). 



An effective Hamiltonian that accurately describes processes which happen on a time 
scale T can now be found by performing the time integrations and identifying the 
dominant terms. For our derivation, we assume that at time t, the excited states of 
the atoms are not occupied and that no photons are present, 

(e|^,t) = and ak\^,t) = 0. (A.4) 

Furthermore, we assume for the parameters in H{t) 



2A. 



< 1 



2A,(A, - 51) 



for X = a,b 



^ 1 for x,y = a,b 



(A.5) 



(A.6) 



We keep terms up to n = 3 on the right hand side of equation (|A.3|) and, by virtue 

IS 



of equation (A.5), neglect all oscillating terms to arrive at 

4 



|^,t+T) 




/i! 



J2^-^^^-^H.T]\^,t),{A.7) 
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N 



Hi 



4A 



Ho 



N 



4A? 



2Afe ' 4(A,-Afc 



+ 



4A, 



2 / 1^ 12 



+ 



j,k I 



4A2 



+ 



2Aa ' 4(A,-AJ ' ^ \6l 



+ 



\9b. 



j,k\ 



Aa 51 - A„ 



ZLA A, 



9b,j,k9a,Lk 



|6j-)(aj-| ® |6;)(az| +H. 



l^bP 9l,j,k9a,l,k 



\bj){aj\ ® |aj)(6/| +H.C. 
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I (^j ) ('^i I 



Provided we can choose the time scale T such that ||iJo^|| ^ Ij ll-f^i^ll ^ 1 and 
||i727'|| ^ 1 but A^T ^ 1 and (5^T 1 (x = a, 6), we can write an effective 
Schrodinger equation, 



\<iI,t + T)-\^,t) 
T 



-I iHo + Hi + H2)\^,t). 



On time scales T, the dynamics of our system is thus accurately described by the 
effective Hamiltonian H = Hq + Hi + H2, which is identical to the effective spin 



Hamitonian of equation (30) up to an irrelevant global constant. 



Appendix B. Explicit expressions for the spin parameters in the N = 2 



case 



B\ 



N=2 



6 
2 



' 8AI 
\9bm 



A 



2A, 4(A, - A5) 



(B.l) 



+ 



(56 

8AI , 



A„)2 



w 

2A, 4(Afe - A, 
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\ga\^{6a-Aa) 



\gb\^{6b-Aa) 



l\N=2 



\9a\'J 



+ 



I a 



4A2 {5a - A,)2 - J2 4A2 (5, - A,)2 - J2 ' 

Here, 5a = oJe — coc and 6b = uje — (i^b — 5) — wc- 



(B.2) 
(B.3) 



Appendix C. Master equation for the dynamics of two three-level atoms 

The three levels of ^"^Rb depicted in figure [6] are two of the 5s^S'i/2 ground-states, 
\a) = \F = l,mp = 1), \b) = \F = 2,mp = 1), and the ^p^Pz/2 excited state 
|e) = |F' = 2, mi?/ = 2). In order to compute the dynamical evolution of this system 
in the presence of light, one needs to include spontaneous emission on the vr 
transition to the 5s^S'i/2 ground state \x) = \F = 2,mF = 2), which is outside 
the three-level system under consideration. The dynamics is thus described by the 
master equation. 



TIT 1 , ^ 

7 
6 

+ 3 Z^^2a^- paj -<yj p- pcTj 

j 



j 

+ fi Z^^2a^- pa^ -<^j P- P(^j ) 



+ K:^^(2ajpaj — aja^p — pa^aj), 
j 

where aj^ = \aj){f3j\ and the atomic transition rates are weighted according to the 
squares of the respective dipole matrix elements . 
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